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We discuss the enumeration of planar graphs using bijections with suitably decorated trees, 
which allow for keeping track of the geodesic distances between faces of the graph. The 
corresponding generating functions obey non-linear recursion relations on the geodesic 
distance. These are solved by use of stationary multi-soliton tau-functions of suitable 
reductions of the KP hierarchy. We obtain a unified formulation of the (multi-) critical 
continuum limit describing large graphs with marked points at large geodesic distances, 
and obtain integrable differential equations for the corresponding scaling functions. This 
provides a continuum formulation of two-dimensional quantum gravity, in terms of the 
geodesic distance. 
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1. Introduction 

The work presented in this note was initially motivated by the need to better un- 
derstand on a combinatorial level the various results obtained via matrix models on the 
enumeration of graphs with fixed topology jl] (see also @ and [|] and references therein) . 
The early work on this subject dates back to a combinatorist, W. Tutte 0, who managed 
to enumerate many of the planar versions of these using recursion relations in the spirit of 
what we call today "loop equations" for matrix models. Such an approach, though com- 
binatorial, failed to really explain the simplicity of the algebraic equations determining 
the generating functions for planar graphs. Another motivation comes from the physics of 
two-dimensional quantum gravity. At the discrete level, coupling matter to gravity simply 
amounts to define a statistical model (typically with local Boltzmann weights) on a fluc- 
tuating base space, in the form of random discretized surfaces. The continuum version of 
this involves field-theoretical descriptions of random surfaces with critical matter |J, via 
the coupling of conformal field theories (matter) to the Liouville field theory (metrics of 
the underlying space). 

The interpretation of the planar graph results remained elusive until the groundbreak- 
ing work of G. Schaeffer ||, who finally gave a beautifully simple combinatorial explanation 
for these algebraic equations, in terms of decorated trees. The idea was simply to estab- 
lish bijections between classes of planar graphs and suitably decorated trees, then easily 
enumerated via algebraic relations obeyed by their generating functions. This technique 
proved quite general, and was extended to many classes of planar graphs, including graphs 
of arbitrary valence , special classes of bipartite graphs called constellations M M and 



other classes of bipartite graphs |T0|], including the particular cases of hard objects on 
planar graphs and of the Ising model on planar graphs ||12|1 . 

The great advantage of this bijective enumeration is that it allows for keeping track of 
some details of the graphs in the language of trees. An important example of this concerns 
the geodesic distance between say the vertices of random quadrangulations, namely planar 
graphs with only tetravalent faces. In ||13|| , it was shown that the geodesic distance of all 
vertices from an origin vertex of the graph may translate into integer vertex labels in some 
corresponding trees, themselves realizing a discrete version of the Brownian snake, further 
studied and extended in |14| in the language of spatial branching processes. In ||15|1 , it 
was shown that the generating functions for planar graphs with two external legs obeyed 
non-linear recursion relations on the maximal geodesic distance between the legs, and it is 
the purpose of this note to clarify and extend the results of this paper. 
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Note that no continuum (field-theoretical) treatment of two-dimensional quantum 
gravity in terms of geodesic distance is available yet, only some partial results were obtained 



using a transfer matrix formalism [16] [17] [18] leading to some conjectural continuum 
expression for a scaling function of the geodesic distance in the case of pure gravity without 
matter. The present work provides an alternative solution and gives access to a host of 
scaling functions for various (multi-) critical models of matter systems coupled to two- 
dimensional quantum gravity. 

The paper is organized as follows. In Sects. 2,3 and 4, we recall some known facts on 
the bijective enumeration of planar graphs with respectively even valences, arbitrary (even 
and odd) valences, and bicolored vertices. This relies on an iterative cutting procedure 
which, starting from a planar graph, produces a decorated rooted tree. In turn the tree may 
be closed back in a unique way into a planar graph, and we use this bijection to recover 
the algebraic relations satisfied by the various counting functions involved. Sect. 5 is 
devoted to the introduction of the geodesic distance in the various enumeration problems at 
hand: we present a unified picture involving formal operators Q generating the descendent 
trees around a vertex, and allowing for keeping track of the geodesic distance from the 
root to the external face of the planar graph to which the tree closes back. The main 
results are recursion relations on the geodesic distance satisfied by the generating functions 
for planar graphs with legs. In Sect. 6, we solve exactly a number of these recursion 
relations, by expanding the generating functions at large maximal geodesic distance n, 
and resumming the resulting series. These display a remarkable "integrable" structure 
in that they generically involve tau- functions of the KP hierarchy. Sect. 7 is devoted to 
the critical continuum limit of the problem, in which we consider large planar graphs and 
large geodesic distances. The generating functions calculated in Sect. 6 are shown to yield 
universal scaling functions, characteristic of the various (multi-) critical points of random 
surfaces with matter and with marked points at a fixed geodesic distance. We propose a 
generalization of our results based on differential equations obeyed by the scaling functions, 
and illustrate it in the case of the Ising model on random surfaces. Finally, we gather a 
few concluding remarks in Sect. 8, where in particular we discuss the integrable structure 
of our recursion relations, and of their continuum counterparts. 



2. Planar graphs and trees I: the case of even valences 

In this section, we recall some results on the enumeration of planar graphs of even 
valence with two extra "legs". This is done via a bijection between planar graphs and 
decorated trees, easily enumerated. 
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(a) 



(b) 





(c) 



(d) 



Fig.l: The bijection between planar tetravalent graphs with two legs (a) 
and rooted ternary blossom-trees (d) is obtained via the following cutting 
procedure. We first visit all edges of the graph in counterclockwise direction 
(b) and cut them iff the resulting graph is still connected: we have indicated 
the succession of cut edges by their number in the order of visit, from 1 to 6 
here. The cut edges are then replaced by pairs of white and black leaves (c) . 
Moreover the incoming leg is replaced by a white leaf and the outcoming one 
by a root, finally yielding a rooted ternary blossom-tree (d). 

2.1. Tetravalent graphs 



We wish to compute the generating function R = R{g) for planar tetravalent graphs 
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with two distinguished (say in- and out-coming) univalent vertices, and with a weight g 
per tetravalent vertex. These will be referred to as "two-leg diagrams" in the following, 
the legs denoting simply the edges connecting the two univalent vertices to the graph. For 
definiteness, we will always represent these planar graphs with the incoming leg adjacent 
to the external face. Note that the outcoming one need not be adjacent to the same face, 
as in the example of Fig.l (a). 

The computation of R relies on the following bijection, illustrated in Fig.l, between 
two-leg diagrams and so-called rooted blossom-trees @. Starting from a two-leg diagram 
(Fig.l (a)), let us first visit all edges adjacent to the external face, starting from the 
incoming leg and in counterclockwise order. Successively, each visited edge is cut iff the 
cut diagram remains connected (Fig.l (b)). The cut ends of the edge are then decorated 
respectively with a black and a white leaf. Once all edges adjacent to the external face 
have been visited, we repeat the procedure with the newly cut diagram. The process 
ends when all faces of the original diagram have been merged with the external one. The 
resulting graph is nothing but a planar tree (it has only one face), with black and white 
leaves (Fig.l (c)). Note that by construction there is exactly one black leaf connected to 
each internal vertex of the tree. Finally, we replace the incoming vertex by a white leaf 
and the outcoming one by a root, so that there is exactly one more white leaf than black 
ones (Fig.l (d)). We define rooted blossom-trees as planar rooted ternary trees with black 
and white leaves, and such that there is exactly one black leaf at each internal vertex. Our 
cutting procedure has produced a rooted blossom-tree out of any two-leg diagram. The 
process however is readily seen to be invertible as there is a unique way of re-connecting 
the black- white leaf pairs into edges, by connecting each black leaf to the first available 
white leaf in counterclockwise order around the tree (the dashed lines of Fig.l (c)). This 
establishes the desired bijection between the two objects. 

Couting rooted blossom-trees is now an easy task, performed for instance by inspection 
of all possible environments of the vertex connected to the root. This leads straightfor- 
wardly to the relation 




R = 1 + 3gR 2 

(2.1) 
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where the first term 1 counts the possibility that the root be directly connected to a leaf, 
and the second term accounts for the three possible positions for the black leaf around 
this vertex, itself receiving the weight g, while the two other descendents of the vertex are 
themselves rooted blossom-trees. From its very definition as counting function, R admits 
a power series expansion in g, with R = 1 + 0(g). This fixes it uniquely to be 



The series for R has a finite convergence radius g c = 1/12. When g approaches g c (critical 
limit), the contribution of large graphs becomes dominant, and we learn that the number 
of graphs with N vertices behaves as g~ N /iV 3 / 2 for large N. 

Note that R(g) = C(3g) where C denotes the generating function for Catalan numbers 
cn = ( 2 ^)/(N + 1), which count among other things the rooted planar binary trees with 
N inner vertices, with the convention that cq = 1. The number of rooted blossom trees 
with N inner vertices is obtained by considering rooted planar binary trees with all leaves 
white and by decorating each vertex with a black leaf: it reads therefore 3 N cjy as there 
are three choices for the position of the black leaf at each inner vertex. 

2.2. General case of graphs with even valences 

The bijection of previous section may be extended to include two-leg-diagrams of 
graphs with arbitrary even valences v = 4, 6, 8... We repeat the exact same cutting proce- 
dure and end up with some generalized rooted blossom-trees, such that each inner vertex 
say of valence v = 2k has exactly k — 1 black leaves attached to it. The corresponding gen- 
erating function R = R(gi, g&, g%, ...) with say weights gik per 2/c-valent vertex, k = 2, 3, 4... 
obeys the following relation 

i?=i+£^( 2 \ _i y (2.3) 

k>2 ^ ' 

obtained again by inspecting all possible environments of the vertex attached to the root. 
The combinatorial factor ( 2k ^ 1 ) accounts for the number of choices for the positions of the 
k — 1 black leaves around the vertex, while the remaining k descendents are themselves 
rooted blossom-trees. 

Again, R is the unique solution to ( |2.3| ) that admits a power series expansion in the 
g 2 iS, with R = 1 + 0(g 2 i) for all i > 2. 
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3. Planar graphs and trees II: arbitrary valences 

This section extends the results of the previous one to graphs with both even and 
odd valences. The first consequence of allowing for odd valences is the existence of one-leg 
diagrams with only one (outcoming) external leg. These will be represented in the plane 
like two-leg diagrams, but the leg need not be adjacent to the external face. 



3.1. Trivalent graphs 

Let S = S(g) and R = R(g) denote the generating functions for respectively one- and 
two-leg diagrams of trivalent planar graphs, with a weight g per trivalent vertex. Applying 
the cutting procedure of previous sections to one and two-leg diagrams, we end up with two 
types of rooted blossom trees which we call S-trees and R-trees respectively. Note that the 
unique (outcoming) leg of the one-leg diagrams is replaced by a root in the corresponding 
blossom-tree. To characterize S- and R-trees, let us introduce the charge q of a tree as its 
number of white leaves minus that of black ones. For instance, the blossom-trees of Sects. 
2.1 and 2.2 above have all charge q = 1, the same holds for the present R-trees, while 
the S-trees are neutral, with charge q = 0. Now S-trees and R-trees are characterized 
among rooted binary blossom-trees as having only descendent subtrees (not reduced to 
black leaves) of charge or 1, while their total charge is and 1 respectively (this is easily 
done by following the effect of the cutting procedure on the original graph, see @ for 
details). The cutting procedure establishes a bijection between one- and two-leg diagrams 
and S- and R-trees respectively. 

The latter are easily counted, again by inspection of the local environment of the 
vertex attached to the root. We get the coupled relations: 
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S = 2gR + gS< 

44 
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+ 9 



A + 9 A 

JOGS) (s)Cr 



R=l + 2gRS 



(3.1) 
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respectively displaying contributions from a vertex with total charge (with one black leaf 
with q = — 1 and one descendent R-tree with q = 1, and two possible positions for the 
black leaf, or with two descendent S-trees), and from a vertex with total charge 1 (with 
one descendent S-tree with q = and one descendent R-tree with q = 1, and two possible 
relative positions for these). 

The generating functions i?, S are uniquely determined by the relations ( |3.1|) and the 
fact that they admit power series expansions R = 1 + O(g), S = O(g). 

3.2. General case 

The trivalent case is easily extended to the case of arbitrary (even or odd) valences 
weighted by (73, f?4, (75, ••• per tri-, tetra-, penta-,... valent vertex, in which the very same 
cutting procedure now leads to generalized rooted blossom S- and R-trees, now blossom 
trees of arbitrary valences v = 3, 4,5, ... again further characterized by the fact that all 
their descendent subtrees not reduced to a black leaf have charge or 1, and by their total 
charge and 1 respectively. This allows to count them straightforwardly, with the coupled 
relations: 

rfe— li 



fc>3 j=0 



3 J \ 3 



(3.2) 



k>3 j=0 

where the combinatorial factors account for the possible ways of positioning j black leaves, 
j or j+1 descendent R-subtrees and the remaining S-subtrees on the vertex attached to the 
root. Again, eqs.(|3?2|) determine completely R and S with R = 1 + 0(gi) and S = O(gi) 
for all i > 3. 

For illustration, in the case of tri/tetravalent graphs, where only #3, #4 are non-zero, 
we have the equations 

S = g 3 (2R + S 2 ) + g 4 (QRS + S 3 ) 

(3.3) 

R = l + 2g 3 RS + 3g 4 R(R + S 2 ) 



4. Planar graphs and trees III: bipartite graphs 

We now turn to the slightly more involved case of bipartite (i.e. vertex-bicolored, say 
black and white) graphs. In the language of matrix models, these correspond to the case 
of two coupled matrices. 
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4-1- p-valent case 

Let us consider two-leg diagrams of vertex-bicolored p-valent planar graphs, and their 
generating functions with a weight g (resp. g) per p-valent black (resp. white) vertex. 
We must also indicate the color of the vertices to which the in and out-coming legs are 
attached, and this leads to a priori distinct generating functions. For simplicity, we restrict 
ourselves to only diagrams with incoming (resp. outcoming) leg attached to a white (resp. 
black) vertex with generating function R = R om (g,g), and also include the single graph 
whose incoming and outcoming legs are directly attached to one-another, without vertex, 
contributing 1 to R. 




(a) (b) (c) 

Fig. 2: A sample bipartite p-valent two-leg diagram with p = 4 (a) is iter- 
atively cut into a rooted bipartite blossom-tree (c) by the usual procedure, 
with the restriction that only edges originating from a black vertex may be 
cut. The edges to be cut are indicated in (b) with their order of visit, coun- 
terclockwise around the graph. 

The above cutting procedure is now slightly modified, with the additional constraint 
that an edge may be cut only if it moreover originates from a black vertex (see Fig. 2 for an 
illustration in the case p = 4). This leads to a new kind of blossom-trees (Fig. 2 (c)), with 
bicolored vertices, and such that black or white leaves may only be connected to vertices 
of the same color, while each black vertex has exactly p — 2 black leaves attached to it. 
The tree still has the property of having one more white leaf than black ones (i.e. a total 
charge of q = 1), as the incoming leg is replaced by a white leaf (which is compatible with 
the above rule, as the incoming leg is attached to a p-valent white vertex). Moreover, the 
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root (former outcoming vertex) is attached to a black vertex. The generating function R 
obeys the relations: 

R = l + (p-l)gX 

(4.1) 

X = ~gR^ 

The first line is obtained by inspection of all possible environments of the root: (i) it may 
simply have one white leaf attached to it or (ii) it may have a black vertex attached to it, 
itself with p — 2 black leaves and a descendent rooted blossom-tree of total charge + (p — 1) 
with a black root, and with generating function X. The second line expresses these latter 
trees according to the environment of the white vertex attached to the root, having p — 1 
descendent rooted blossom-trees of charge +1, all generated by R. We finally get 

R=l + (p-l)g~gR p - 1 (4.2) 

Note that the generating function C p (x) for rooted p-valent planar trees with a weight x per 
vertex satisfies C p (x) = 1 + xC p (x) p ~ 1 and C p (x) = 1 + 0(x). The corresponding number 
of trees with iV vertices reads cff = 1+( - p 1 _- L ^ Af (^ P ~^ N ) ■ These numbers are also known as 
the Fuss-Catalan numbers, and reduce to the ordinary Catalan numbers for p = 3. Finally 
the number of rooted blossom-trees with iV vertices is simply R\ g NgN = (p— 1) N C^ ■ For 



p 



large N, it behaves as g c N /iV 3//2 , where g c = (p — 2) p 2 / (p — 1 
4-2. p- constellations 

The p-constellations, introduced in || , are vertex-bicolored (black and white) graphs 
such that say black vertices have all valence p, while white ones may have valences arbitrary 
multiples of p. We again consider two-leg diagrams of such planar graphs with say incoming 
leg connected to a white vertex and outcoming leg connected to a black vertex (see Fig. 3 
(a) for an illustration with p = 3), or both legs being directly connected without vertex, 
and count them with a weight g per (p-valent) black vertex and weights g m for mp-valent 
white vertices, m = 1,2,3... Let R = R°'(g;{gi}) denote the corresponding generating 
function. The cutting procedure, illustrated in Fig. 3, remains the same as in the previous 
section, and leaves us with rooted blossom-trees with bicolored vertices of total charge +1, 
such that leaves may only be connected to vertices of their own color, and that descendent 
subtrees may be of two types: if their first vertex is black, they have total charge 1; if its is 
white, they have total charge p — 1, the descendents of the root vertex of the latter being 
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1 



(a) 



(b) 



(c) 



Fig.3: The cutting procedure is applied to a two-leg diagram of a 3- 
constellation (a), with incoming (resp. outcoming) leg attached to a white 
(resp. black) vertex. The edges are visited in counterclockwise order around 
the graph (b), and cut iff (i) this leaves the resulting graph connected (ii) the 
cut edge originates from a black vertex. We have indicated the chronological 
order of the cut edges, from 1 to 6. We finally replace each cut edge by a pair 
of black/white leaves, and the incoming (resp. outcoming) leg by a white leaf 
(resp. root). 



themselves only blossom trees of charge 1 or bunches of p — 1 black leaves attached to a 
black vertex. 

By inspection of all possible local environments of the vertex attached to the root, we 
may derive the following relations for R: 



where R, X , Y generate rooted blossom-trees of respective charges 1, p — 1, 1 — p, starting 
respectively with a black, white, black vertex. Y = g is due to the fact that the only tree 
contributing to Y has a black vertex and p — 1 black leaves attached to it. In the case 



R 



l + (p 




(4.3) 



Y 
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p = 3, this reads 



* + 9 A + 9 A 



® 



R = l + 2gX 



® "J\* 



X = ~ gi R 2 + 5g 2 YR 4 + 



® " A 

Y = g 

Eliminating Y and X from ( f4.3| ), we arrive at the single algebraic equation 

a = i + (p - 1) E ^ ( mp z% m - l R {p - 1)m 

m>l ^ m ' 



(4.4) 



(4.5) 



R is the unique solution to this equation such that R = 1 + 0(g, ^). 

Note that in the case p — 2 of 2-constellations, we recover the even-valent graph result 
fl2.3p, upon taking ^2fe = ^fc for > 2, while ^i = and g = 1- Indeed, in 2-constellations, 
the (2-valent) black vertices may be viewed as decorations of the edges of an arbitrary 
graph with only white even-valent vertices. Setting g = 1 precisely allows to forget about 
these decorations, while g± = simply eliminates 2-valent white vertices. 

4-3. Planar bipartite graphs and the Ising model 

Constellations are easily tractable objects, essentially due to the triviality of one type 
of vertex (black here), whose valence remains fixed. More generally, we would like to 
consider in all generality vertex-bicolored graphs of arbitrary even valences for vertices of 
both colors. For the sake of simplicity, we will restrict ourselves to two- leg diagrams in 
which both legs are attached to white vertices. Let us denote by R = R°°({g2i}; {<72i}) 
the generating function for two-leg diagrams of planar bipartite graphs with weights gn 
(resp. #2i) P er 2i-valent black (resp. white) vertex, % = 1,2,..., and such that both legs 
are attached to white vertices. 
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Applying the now usual cutting procedure to two-leg diagrams, we end up with 
blossom-trees of total charge +1 with bicolored vertices of arbitrary even valences. Their 
characterization however is quite delicate, as it involves describing all their possible de- 
scendent subtrees. These come in two forms according to the color of their vertex attached 
to the root: if the latter is white, the possible descendent subtrees are either reduced to 
a black leaf (of charge —1) or rooted vertex-bicolored blossom trees of total charges 1, 3, 
5, ... ; if it is black, the possible descendent subtrees are rooted vertex-bicolored blossom 
trees of charges 1, —1, —3, —5, ... Let us introduce the generating functions Ri for vertex- 
bicolored blossom-trees of total charge 2i — 1, % = 1, 2, 3, ... and whose root is attached to a 
white vertex, together with the generating functions Xi for vertex-bicolored blossom-trees 
of total charge 1 — 2i, i = 1,2,3,... and whose root is attached to a black vertex, while 
V = X generates vertex-bicolored blossom-trees of total charge 1 and whose root is at- 
tached to a black vertex, or the tree made of a single leaf attached to the root, without 
any vertex (contributing 1 to V). We also introduce the generating function R — 1 for 
the tree made of a single black leaf attached to the root, without any vertex. We now 
simply have to enumerate all possible environments of the vertex attached to the root of 
each of these trees, according to the type of its attached descendent subtrees. This gives 
the following system 



V = l + ^2g 2 k ^2 R h R h- R hk-i 

k>l Jl.J2.--.J2fc-l>° 

Sj;=fc 

X m = ^2g 2 k ^2 R ji R 32--- R 32k-n m = 1)2,3,. 

fc>l Jl.J2.--.J2fc-l>° 
£j£ = k — m 

R m = ^2g 2k ^ X h X h--- X hk-n 771 = 1)2,3, 



(4.6) 



k>l Jl.J2»---»J2fe-l^° 
Sjj = k — m 

The desired generating function R = R\ is the unique solution to this system where all 
.R's, X's and V admit power series expansions of the g's and (?'s. 

The case of the Ising model on planar tetravalent graphs may be viewed as a particular 
case of the above, in which only g 2 , (74, g 2 , cja are non-zero. To see this, recall that the Ising 
model on tetravalent planar graphs is defined by say coloring the vertices of an arbitrary 
tetravalent planar graph in black or white (colors stand here for the spin up or down), and 
counting the configurations with different "nearest neighbor interaction" weights for edges 
connected to vertices of the same color (weight e K ) or of different colors (weight 1), while 
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black (resp. white) vertices are counted with a weight ge H (resp. ge~ H ). Here K and 
H are respectively the spin coupling and the external magnetic field of the Ising model. 
To make the contact with our model, we just have to resum all possible configurations 
obtained by adding arbitrary numbers of 2-valent black and white vertices on the edges 
of Ising configurations, in such a way that bicoloration is restored. This entails adding 
any chain of black, white, black, white 2-valent vertices between any white and black 
tetravalent vertices connected by an edge, or any chain of alternating white, black, 
white 2-valent vertices between any two black tetravalent vertices connected by an edge or 
else any chain of black, white, black 2-valent vertices between any two white tetravalent 
vertices connected by an edge. Doing the resummations within the configurations of our 
bicolored graphs produces an effective edge interaction weight w a b according to the colors 
a, b of the adjacent vertices: 

92 

w„ = - — = •— O— • + •— O— •— O— • + ... 

1 - 9292 
92 

W QO = — = o— • — o + o— • — o— • — o + ... 

1 - 9292 

1 

w.o = w , = — = O— • + O— • — O— • + ... 

1 - 9292 

(4.7) 

hence to identify our model with the Ising one, we must take gi = gi = e K , while g^ = 
(1 — e 2K ) 2 ge H and g± = (1 — e 2K ) 2 ge~ H , and the external legs must receive the extra 
weights 1 / y/l — e 2K each. 

Restricting to the symmetric case gi = g<i = c and g± = g^ = g, the above equations 
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simply read 



® 




/// / /// / 
8 




V = 1 + cR 1 + 3gRl + 3gR 2 




X 1 = c + 3gR! 

/// / /// / 
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x 2 = g 




R 1=C V + 3gV 2 X 1 




R2 = gV- 



(4.8) 



Eliminating R 2 and Xi, we get 

R(c + 3gR) 2 (l - (c + 3gR) 2 ) = (c + 3gR) 3 + 3g 2 R 3 (4.9) 



V = 



with 



c + 3gR 

Here, R = R± is the generating function for two-leg diagrams of planar tetravalent graphs 
with Ising (black or white) spins decorating their vertices, with interaction weights w„ = 
Woo = c, w, = w * = 1 and a weight g/(l — c 2 ) 2 per tetravalent vertex, and such that the 



two legs are attached to univalent black vertices, themselves weighted by 1/yl — c 2 . 



5. Geodesic distances 

In the previous sections, we have enumerated various one- and two-leg diagrams by 
establishing bijections with suitable classes of blossom-trees. Note that in the plane rep- 
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resentation we have chosen, the face Fi adjacent to the unique leg for one-leg diagrams is 
not necessarily the external face F . Accordingly in the case of two- leg diagrams, the face 
Fi adjacent to the out-coming leg need not be the external face F , itself adjacent to the 
in-coming leg. In this section, we show how to keep track of the geodesic distance between 
F and F±, namely the smallest number of edges to be crossed in a path from F to F\. 
By a slight abuse of language, this distance will also be referred to as the distance between 
the legs. 

5.1. Keeping track of the geodesic distances 

The main feature of the previous sections is a sort of unified formulation of planar 
graphs in the language of blossom-trees. Note that going back from blossom-trees to graphs 
is a straightforward step, as there is a unique way of reconnecting the black and white leaves 
into edges: this is done by simply connecting each black leaf to the first available white leaf 
in counterclockwise direction. This process leaves us in the case of trees of charge 1 with 
exactly one unmatched white leaf, which is taken as outcoming leg, while the root is the 
incoming one. For trees of charge 0, all pairs are exhausted and only the root remains as 
uniqe leg. In both cases, we note that the geodesic distance between the faces Fq (external) 
and F\ (adjacent to the former root) is simply given by the number of black-white edge 
pairs that separate the root from the external face after recombination. Keeping track of 
this geodesic distance simply amounts to keeping track of the black leaves "in excess" that 
require encompassing the root to be connected to their white alter ego. 

This is done in all generality by attaching to each blossom-tree a "contour walk", 
namely a walk on the relative integer line with steps ±1, obtained as follows (see Fig. 4 for 
an example). One starts from the root of the blossom-tree and visits in clockwise direction 
all leaves around the tree. Starting from the coordinate 0, we make a step +1 (resp. —1) 
for each encountered black (resp. white) leaf. In a blossom-tree of charge k, such a walk 
will end up at coordinate —k. Now the number of excess black leaves responsible for the 
geodesic distance between F\ and Fq is simply the maximum coordinate reached by the 
contour walk. 

In view of this result, and of the form of all relations determining the blossom-trees of 
the previous sections, it is natural to introduce by analogy with the generating functions 
say X for some particular type of blossom-trees the generating function X n for the same 
blossom-trees with a geodesic distance of at most n between Fi and F . To obtain from 
X n the generating function for blossom-trees with geodesic distance equal to n betwen Fq 
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(c) 



Fig.4: The contour walk of a rooted blossom-tree. Visiting the tree (a) in 
clockwise direction starting from the root, one keeps a record of the type 
of leaves encountered in the form of a walk on the integer line (b), starting 
at the origin, and with steps up (for a black leaf) and down (for a white 
one). The maximum reached by the walk is nothing but the geodesic distance 
separating the root from the external face in the recombined graph. This 
distance is 3 in the present case, as one readily checks by closing the tree 
back into a planar graph (with tetra/hexavalent vertices here). Concentrating 
on the environment of the vertex attached to the root, we see that each 
descendent subtree corresponds to a portion of the walk (c), with a certain 
relative maximum. Expressing the global maximum of the contour walk in 
terms of the relative maxima of its portions allows for writing a recursion 
relation for the generating function for rooted blossom-trees whose root is at 
a maximum distance n from the external face of the recombined graph. 

and Fi, we simply have to take the difference X n — X n _±. Keeping track of n then boils 
down to expressing the maximum of the contour walk of a tree in terms of those of the 
individual contour walks of the blossom-trees descending from the vertex attached to the 
root, following the same inspection procedure as before. This is done case by case in the 
following sections. 

5.2. Even valences 

Let us start with the tetravalent case. We find that eq. 

into 

R n = 1 + gR n (R n -i + R n + Rn+i) (5-1) 



2.1) must be transformed 
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where R n denotes the generating function for tetravalent planar graphs with two legs at 
distance at most n, and weight g per vertex. The three terms on the r.h.s. correspond 
to respectively the black leaf on the right, in the middle or on the left of the two other 
descendents of the vertex attached to the root (see the picture of eq . (|2 ■ 1| ) ) ■ It is clear 
that the presence of the black leaf acts as a shift by —1 on the local distance n while R n 
accompanies a shift by +1, when going clockwise around the vertex. 

To have a compact notation for the general result, let us introduce a formal orthonor- 
mal basis \n), nfZ, with (m\n) = 6 m)n , and an operator a acting as a shift a\n) = \n + l), 
and its formal inverse a -1 such that <j~ l \n) = \n — 1). We introduce the operator 

Q = a + a ~ 1 r (5.2) 

where f simply acts diagonally as f\n) = R n \n). Note that the shift a may be represented 
by a black leaf, while its inverse always accompanies an R. Then eq. ( p.l| ) takes the form 

l = (n-l\(Q-gQ 3 )\n) (5.3) 

More generally, in the case of arbitrary even valences, we have to write 

l = {n-l\{Q-Y,92kQ 2k - l )\n) (5.4) 

k>2 

For instance in the case of tetra- and hexa-valent graphs, eq. (|5.3| ) reads explicitly 

Rn = 1 + g^RniRn+l + Rn + -Rn-l) + 9&Rn (Rn+lRn+2 + Rn+lRn-1 

(5.5) 

+ Rn-lRn-2 + R n +l + R-n-l + -Rn(2-R n +1 + Rn + 2-R n _l)) 

A remark is in order. The equations ( [5.1| )-( |5T5] ) are valid only for n > 0, provided we 
use R-k — 0, k = 1, 2, ... wherever they occur in the r.h.s. With these boundary conditions, 
and the general fact that R n possess a power series expansion in g, with R n = 1 + O(g) 
for n > 0, all i? n 's are then uniquely determined by ( |5.4[ ) order by order in g. 

5.3. Arbitrary valences 

Let S n (resp. R n ) denote the generating function for one- (resp. two-) leg diagrams 
of planar graphs with arbitrary valences with weights gi per i-valent vertex, and such that 
F (the external face) and Fi (the face adjacent to the unique (resp. outcoming) leg) are 
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distant by at most n. For trivalent graphs, we find that eqs. ( |3.1|) must be transformed 
into 

S n = g(R n + Rn-i) + gS n 

(5.6) 

R n — 1 + gR n (S n+ i + S n ) 

where in addition to the situation of previous section we simply note that S"s don't affect 
the distance counting (no shift). This suggests to introduce in the general case of arbitrary 
valences the operator 

Q = a + a~ 1 sa + a~ 1 f (5.7) 

where s acts diagonally on the basis \n) as s\n) = S n \n), and in terms of which we simply 
have to write 



= (n\(Q-J2^Q l - 1 )\n) 

i>3 

1 = (n-l\(Q-J2g i Q i - 1 )\n) 



(5.8) 



i>3 

For tri- and tetra-valent graphs, this reads 

Sn = g%(Rn + Rn-i + S n ) + g4(R n (S n+ i + 2S n ) + R n -i(S n -i + 2S n ) + 

Rn = 1 + g?,Rn{S n + SVi+l) + g4Rn(Sn + S n S n +l + S n+ i + Rn+1 + Rn + Rn-l) 



(5.9) 



Note that when = 0, we find the solution S n = identically, and eq. ( |5.9|) reduces 
to the tetravalent case (|5.1| ). More generally, imposing that all odd g's vanish leads to the 
solution S n = (as there are no one- leg diagrams with only even valences), and we recover 
the even-valent case of previous section. 

5-4- Constellations 

Let R n denote the generating function for two-leg diagrams of p-constellations with a 
weight g per white (p-valent) vertex, and weights cji per black pz-valent vertex, i = 1,2,..., 
whose incoming (resp. outcoming) leg is attached to a black (resp. white) vertex (or both 
are connected without vertex and the corresponding unique graph contributes 1 to R n for 
all n) and such that the geodesic distance between the two legs is at most n. Here the 
notion of geodesic distance is defined according to the rules used in the cutting procedure of 
Sect. 5. 3, namely the geodesic distance between Fq and F\ is the minimal number of edges 
to be crossed in a path going from F to F\ , and such that at each edge-crossing the white 
vertex is always on the right. Following the same reasoning as in the previous sections, we 
are now led to the introduction of two operators Qi and Q2 which generate, upon taking 
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powers, the successive decorations of the vertex attached to the root, respectively in the 
case of a black and white vertex. We have 

Q x = a + a~ 1 xa 2 ~ p 

(5.10) 

Q 2 = a 1 f + a p 2 ya 

where the shift a represents a single black leaf, while x represents rooted blossom-trees of 
charge p — 1 whose first vertex is white. Again, x, y 1 r act diagonally on the basis \n) with 
eigenvalues X n ,Y n ,R n respectively. The equations (J4.3|) now become 



l = (n-l\{Q 2 -gQ p 1 - 1 )\n) 



0=(n + p-l\(Q 2 -gQ p l - 1 )\n) (5n) 
0=(n-p + l\(Q 1 -J2~mQ P 2~ 1 )\n) 



i>i 

In the particular case of 3-constellations, with say only (7,(71,(72 non-zero, these read for 
instance 

R n = l + g{X n + X n - 1 ) 

Y n = g (5.12) 

X n = (}iR n R n+ i + Cj2gRnRn+l{Rn+2,Rn+2 + -Rn+2-Rn+l + Rn+lRn + R n Rn-l) 

while in the case of only (black and white) p-valent vertices of Sect. 4.1, eqs.( |5.1i| ) read 

R n = 1 + g{X n + X n -i + ... + X n - p+ 2) 

(5.13) 

X n — giR n R n +i...R n + p -2 

for the generating function R n for two-leg-diagrams of bipartite p-valent graphs with in- 
coming leg attached to a white vertex and outcoming leg attached to a black one, and such 
that the geodesic distance from the in- to the out-coming leg is at most n. 

5.5. Bipartite even-valent graphs and the Ising model 

In the general case of bipartite even-valent graphs, we are led to the introduction of 
two operators Q%, Q2 with the following structure 



(5-14) 

Q 2 = o-^v + ^o-' 1 



x^a 2k 



k>i 
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where the A k \x^ k \v all act diagonally with eigenvalues R n \x n ,V n . The latter are 
nothing but the generating functions for sets of rooted blossom-trees restricted by n, re- 
spectively starting with a white, black, black vertex, and with charges 2k — 1, 1 — 2k, 1 
respectively. We are actually interested in computing R n = Rn , the generating function 
for two-leg diagrams of bipartite planar graphs with weights g 2 i,g 2 i per 2z-valent black, 
white vertex, such that moreover the two legs are distant by at most n. Again, the distance 
from Fq to F\ is defined as the minimal number of edges to be crossed in a path from Fq to 
such that at each edge crossing the white vertex is always on the right. This definition 
allows to keep track of this distance on the trees themselves, as the number of "excess" 
black leaves which upon recombination with white ones encompass the root of the tree. 
The equations determining R n are simply 

l={n-l\{Q 2 -Y J 92iQl i - l )\n) 
i>i 

0=(n + 2m- 1|(Q 2 - ^^<5f 

i>l 

= (n - 2m + l|(Qi -J^^iQf'^W), 

i>l 

In the abovementioned case of the Ising model with only g 2 = g 2 = c and g& = = g 
non-zero, we must take Qi = a+a -1 ?^ +a~ 3 r^ and Q 2 = a~ 1 v+a~ 1 x^a 2 +a~ 1 x^a 4 7 
and eq.( |5.15| ) reduces to 

V n = l + cR n + gR n (R n+1 + R n + R n -t) + g(R^ + R^ +1 + R%) 

=C + g(R n + Rn-l + Rn-2) 

X^=g (5.16) 
Rn = cV n + gV n (V n+1 X n % + Ulji + Vn^xW) 

Rn 2) = gVnVn-xVn-2 

We first remark that R n = VnX^h by comparing the second and fourth lines of eq.( [5.16| ). 
This is a particular case of a general duality between Qi and Q 2 in the symmetric case 
when $j = gi for all i, where we may write Qi = Q\, where ert = cr~ l v , (AB)^ = B^ A^ for 
all operators A, 5, and f* = f for all diagonal operators. Here this implies 
r^v-V and a~ x x^a A = r^v" 1 *) 3 , i.e. R n = V n X n % and R n 2) = 
Finally eliminating X^ and R n 2 ^ from eq.( |5.16| ), we are left with 

V n (l - g 2 (V n+1 V n+2 + Vn+xVn-x + V n - X V n - 2 )) = 1 + R n (c + g(R n +i + R n + R n -i) 
R n = V n (c + g(R n+ i + R n + R n -i) 

(5.17) 



m = 1,2, ... 



(5.15) 



m = 1,2, 
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An important remark is in order about the generating function R n . Although n has the 
meaning of a maximal geodesic distance between the two legs in the bipartite graph picture, 
it loses somewhat of its meaning in the correspondence with Ising model configurations. 
Indeed, within configurations of the Ising model on tetravalent planar graphs, n is not 
the obvious geodesic distance between the faces adjacent to the in- and out-coming legs, 
as its definition involves first transforming the graph into a bipartite one, and it then 
corresponds to a distance where edge-crossing is permitted only if the white vertex is on 
the right. This restriction is almost irrelevant, as there are in general sufficiently many 
successions of black, white, black... bivalent vertices to allow for crossing edges in both 
directions. One case however is troublesome: when an Ising edge connects a black vertex 
to a white one, in the absence of intermediate bivalent vertices (not necessary here as the 
bicoloration is already ensured) , the edge may only be crossed in one direction, leaving the 
white vertex on the right. This introduces a bias in the notion of distance, having to do 
with the matter configurations on the graph. An analogous situation was encountered in 
TUl [15] in the case of hard dimers on tetravalent planar graphs, namely on configurations 



of tetravalent planar graphs where edges may (or may not) be occupied by dimers which 
repel one-another in such a way that no two adjacent edges can be simultaneously occupied. 
In this case indeed, the notion of geodesic distance is biased by the dimers, in that the 
occupied edges cannot be crossed in paths from Fo to F\. In both Ising and hard-dimer 
cases, the matter interfers with the space, by modifying the rules governing distances. 
Consequently, we think it is interesting to investigate the dependence of the Ising two-leg 
diagrams on this special distance, and it may eventually be that its difference with the 
true geodesic distance becomes irrelevant in large graphs. 



6. Exact solutions 

6.1. Finding exact solutions: a general scheme 

All the equations listed in Sect. 5, despite their diversity, are all basically of the same 
form: (possibly coupled) algebraic recursion relations expressing R n (and the other generat- 
ing functions involved) in terms of a finite number of previous terms _R n _i, _R n _2, R n —k- 
In principle the boundary data needed to entirely determine R n should consist of k con- 
secutive initial values of Rj. It turns out however that we may drastically simplify, namely 
divide by 2 this required number of initial data by simply requiring that lim n ^ 00 R n exists, 
and that it moreover coincides with the generating function R. Indeed, this is nothing but 
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restating the definitions of R n and R, as the latter was first obtained regardless of the 
geodesic distance between legs, while the limit n — > oo of the former amounts to removing 
the geodesic distance constraint in the counting of graphs. In the same fashion, all other 
generating functions involved tend to their obvious limiting values when n — > oo. This 
allows to linearize the various recursion relations at large n, by setting say R n = R — p n , 
and similarly for the other generating functions involved. At first order in p n and its 
other counterparts, we obtain (possibly coupled) linear recursion relations. We immedi- 
ately deduce that p n ~ x n (or a linear combination involving x's of the same modulus) for 
some solution x (with modulus less than 1) to the characteristic equation of the linearized 
recursion relations. 

To completely solve our equations, we start by determining the exact form of the lin- 
earized solution, in general a linear combination p n ~ X/?=i a j( x j) n where Xj, j = 1, 2, k 
denote all the (generically distinct) solutions of the linearized characteristic equation with 
modulus less than 1. In a second step, we obtain order by order in the (xj) n the higher 
order contributions to the true solution, expanded at large n. These take in general the 
form of recursion relations for the coefficients of the multiple expansion in powers of the 
(x-,) n . Solving these recursion relations, and resumming the resulting series allows us to 
finally obtain compact expressions for the exact solutions to the non-linear recursion rela- 
tions at hand. The result still depends on the initial parameters a,j, j = 1, 2, k, which 
then are fixed by requiring that the terms involving the k first R-k, R-k+i, -R-i drop 
out of the recursion relations. 

In the following sections, we simply present the solutions, as we have found them. A 
case by case proof by substitution is left as an exercise to the reader. In many situations, 
the proof boils down to a certain identity between Chebyshev polynomials of the first kind, 
as will be apparent soon. 

6.2. Tetravalent case 

For pedagogical purposes, we detail in this simple case the general scheme presented 
in the previous section. Substituting R n = R — p n into eq. (|5.1|) we get at first order in p n : 




(6.1) 



The linearized characteristic equation therefore reads 




(6.2) 
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For g < g c = 1/12, there is generically a unique solution x = x(g) with modulus less than 
1 to this equation, and we find that at first order p n = ax n + 0{x 2n ). We may now infer 
the general form p n = J2j>i CLjX 71 ^ , Oi = a for the complete solution, where the coefficients 
dj are to be determined order by order in x n . We find explicitly 

k ( x J + i + 1 \ 
a k+ i = [ x k^ , T r ) Wk+i-j (6-3) 



solved recursively as 



1 — x ( ax 



ak=a T^x-{(l-x)(l-x*)) (6 - 4) 
Picking a = x(l — x)(l — x 2 )A, R n is easily resummed into 

R n = R UnUn+S , u n = l-Xx n+1 (6.5) 

Un+lUn+2 

This is the general solution to eq. (|5.1|) , that converges for large n. To see why, it is simplest 
to substitute the form ( |6.5| ) into the initial equation ( |5.1| ) which then boils down to the 
following quartic relation 

u n u n+1 u n+2 u n+3 = ^ul +1 ul +2 + gR(u n _xul +2 u n+3 + u 2 n u 2 n+z + u n u 2 n+1 u n+A ) (6.6) 

Substituting u n = 1 — Xx n+1 into this, we just have to check that the zeros of the l.h.s. 
as a degree 4 polynomial in A = Xx n match those of the r.h.s. as moreover the equation 
reduces to eq. ( |2.1| ) for A = 0. 

Finally, the "integration" constant A is now fixed by further requiring that eq. (|5.1|) 
makes sense at n = 0, in which case the term must drop off the r.h.s. of the recursion 
relation, in other words we have to impose R-% = 0. This simply gives A = 1, and finally 
the exact solution to our combinatorial problem reads 

d ^ ff (l-x"+ 1 )(l-x-+ 4 ) _ U n U n+3 ( 1 \_ x^-x-^ 

n K (i-^+ 2 xi-*-+3) K u n+1 u n+2 > Un - Un [y x+ ^j x \- x -\ 

(6.7) 

where we have identified the Chebyshev polynomials U n of the first kind. 

As a simple application, the formula (|6.7D gives access to the generating function 
for two-leg diagrams whose legs lie in the same face, already identified as that of rooted 
tetravalent planar maps or quadrangulations. We find that 

R Q = R X+ * = R l ~ AgR =R-gR3 (6.8) 

with R as in (|2.2|) . This result was first obtained by Tutte [Q in a completely different, 
though combinatorial, manner. 
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6.3. Even valences 

Let us for definiteness consider the equation ( |5.4|) with only g 4 , g 6 , g2 m +2 non-zero. 
Linearizing again the equation at large n and solving for the leading p n ~ x n , we find that 
x must obey the following characteristic equation: 

'2k + l\ 1 1 - x 2k+1 - 21 

x k-l 1 — x 

(6.9) 



o = Xm (x) = i-jr 92k+2 R k j2( 2k t 1 ) 

I — n 7 — n V / 



k=0 1=0 
m k 



k=l 1=0 ^ ' 



expressed in terms of Chebyshev polynomials of w = y/x + l/y/x. Note that Xm{ x ) * s 
actually a degree m polynomial in x + l/x, with generically m distinct roots with modulus 
less than 1, denoted by x±, x 2 , x m . Repeating the straightforward, though tedious, 
exercise of previous section, we end up with the following exact solution 

(m) (m) 

r>n (m) (m) 

U n+l U n+2 

m I 

«i m) = £(-i)' E n A -C II ^ ( 6 - 10 ) 

1=0 l<mi<...<mi<mi=l ^_i<3_l 

= (x a -x b ) 2 

Ca,b — 



1 X a Xb) 



Remarkably, the structure of «l m ' matches exactly that of an iV-soliton tau-function of 
the KP hierarchy [|19| which reads 



C-4 , 



TV r / 

r=0 ii <...<i r /J=l \/i<f 

_ (Pi-Pj)(Qi ~Qj) 



(6.11) 



' l ' J (Pi ~ Qj)(Qi - Pj) 
where the 77's contain the times' dependence of the KP hierarchy. Our solution ( |6.10|) sim- 
ply amounts to identifying N = m, e Vi = — x™ +m Ai, pi = Xi and % = 1/xi, i = 1, 2, m. 
This surprising relation suggests the existence of an underlying integrable structure for our 
recursion relations, but remains mysterious. 

Again, to further fix the "integration constants" Ai, A2, A m , we simply have to 
express that U-\ = U-2 = ■■■ = U- m = 0, with the result 

A . = Trl^^i z = l,2,..,m (6.12) 

-L -L rn . nf . 

. , . J^l ^ 7 
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and finally the complete solution to (|5.4j) reads 



U n (w 1 ,...,w m )U n+3 (w 1 ,...,w m ) 
R n = Rj- T 7 7 c (6.13) 

u n+ i{wi, w m jty n+ 2(wi, ... ,W m ) 

where 

U n (wi, u> m ) = det [f/ n+2 j-2(^i)]i<i,j< m (6.14) 
in terms of Chebyshev polynomials of the first kind expressed at the values Wi = ^fxl + 



1/y/xi, % = 1, 2, m. The precise proof of these statements can be found in [[15 . 

For illustration, in the tetra/hexavalent case m = 2, we have the characteristic equa- 
tion 



= X 2(x) = 1 - g 4 R ( x + - + 4 J - g Q R 2 ( x 2 + 6x + 16 + - + ) (0. 15) 



and for instance the result ( |6.13| )( |6TT4D reads for n = 0: 



1 + + -L)(l + I2 + J_) 1 - 3g 4 fl - 1096-R 2 V 

where we have reexpressed the symmetric functions of X\ + 1/xi and X2 + 1/^2 in terms 
of the coefficients of Xi{x). The result is nothing but the generating function for rooted 
tetra/hexavalent planar graphs. 

6.4- Trivalent case 

Repeating the exercise of Sect. 6. 2 for the trivalent case of eqs. flSDD , we have found 
the following solution 

Rn = R }^±l, u n = l-Xx n+1 (6.17) 



u l+i 



while 



S n = S — — , t n = gR 2 (l-x)(l-x 2 )Xx n (6.18) 

U n U n +l 

In eqs. (|6.17|) ( |6.18| ), the parameter x is the unique solution to the linearized characteristic 
equation 

l-g 2 R 3 (x + - + 2 J =0 (6.19) 



x 



such that \x\ < 1. 
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Requiring that R-i = 0, fixes A = 1, and we obtain the complete solution 

,(1 - x n+1 )(l - x n+3 ) 



Rn — R- 



S n = S-gR 



(l-x n + 2 ) 2 

(1 — x)(l — x 2 )x 



( 6 - 2 °) 



(1 - x n+1 )(l - x n + 2 ) 



from which we read off the following compact expressions for R and So, respectively 
generating two- and one-leg diagrams of trivalent planar graphs, with respectively the two 
legs in the same face and the leg in the external face: 

x + x + 1 2 4 

-Rn = R T = R — Q R 

x + I + 2 (6.21) 

S = S-gR 2 

In particular, R is the generating function for rooted trivalent planar graphs. 

6.5. Arbitrary valences 

The case of arbitrary valences still awaits a good solution, in the same spirit as the 
case of even valences. It is however possible to find "integrable-like" solutions to the 
corresponding recursion relations (containing only one integration constant) , which for the 
time being are still too restrictive to describe the general case. Indeed, we need a sufficient 
number of integration constants to allow for satisfying all the necessary initial conditions 
of our combinatorial problem. This number is exactly the degree of the characteristic 
equation of the linearized recursions, when expressed as a polynomial in x + 1/x. 

For arbitrary values of #3, g±, g$, ... we have the following "one-x" solutions: 

(1 - Xx n+1 )(l - Xx n+3 ) 

Rn — R — ^ „_i_o\ 2 



Sfi — *S* v Rx 



(1 - Xx n + 2 ) 

(1 - x)' 2 Xx 



(6.22) 



(l-Ao:"+ 1 )(l-A^+ 2 ) 



where x is any solution with modulus less than 1 to the corresponding linearized charac- 
teristic equation. The latter reads for instance in the case of tri/tetravalent graphs where 
only gs and #4 are non-zero: 

g 4 R(x + -+4) + S(2g 3 + 3g 4 S)-l) - R(g 3 + 3g 4 S) 2 (x + - + 2) = (6.23) 
x J x 

26 



and amounts to 

\/%3 + 3g 4 S) {y/i + -L) = l - g 4 R( x + - + 4) - S(2g 3 + 3g 4 S) (6.24) 



as R has a power series expansion of the form R = 1 + 0{g 3l g 4 ). In this particular case, 
it would be desirable to obtain the full solution involving the two roots x\ and x<i of 
( |6.23|) and two integration constants, to be able to solve simultaneously the two boundary 
conditions R-i = and R-iS-i = (to be understood as lim n ^_i R n S n = 0) obtained 
from (|5.9f) at n = 0, and clearly not satisfied by ( |6.22|) . 

6.6. Bipartite p-valent case 

Repeating the usual exercise with eq.( |5.13|) in which we first eliminate X n , we have 
found the following solution, including only one integration constant: 

R n = R UnUn +P+ 1 , Un = i _ Ax n+1 (6.25) 

valid provided x is chosen among the roots of modulus less than 1 of the linearized char- 
acteristic equation: 

1 - gg 1 R p ~ 2 -^(l + x + x 2 + ... + x p ~ 2 ) 2 = (6.26) 

As noticed in the case of planar graphs of arbitrary valences, the degree of this polynomial 
of x + 1/x is however p — 2, hence only for p = 3 is the solution ( |6.25| ) completely general. 
This is the case of bipartite trivalent planar graphs, for which eq.( |5.13j ) reduces to R n = 
1 + ggiR n (R n +i + R n -i). In this case, the solution is further fixed by requiring R-i = 0, 
hence A = 1, and it reads 

with 1— ggiR(x + l/x + 2) = 0, \x\ < 1. R n is the generating function for bipartite trivalent 
graphs with two legs attached to vertices of opposite colors, and geodesic distance between 
those less or equal to n. In particular, 

Ro = H^±4+£+i±l = fii-frftfl+ffji? (628) 

(x + i)(x+J+2) l- 2 99l B 

is the generating function for rooted bipartite trivalent planar graphs (with weights g/gi 
per black/white vertex) 



27 



For p > 4, we must work out the generalizations of ( |6.10|) , which now read as follows. 
Introducing 

p(x) = x(l + x + ... + x p_2 ), q(x) = p(l/x) 

Pi=p(xi), qi=q(xi) / 62g x 
. (Pi ~ Pj)(Qi ~ Qj) 



C i,3 



(Pi - Qj)(Qi-Pj) 

where Xi, i = 1,2, ...,p — 2 denote the generically distinct roots of eq.( |6.26|) with modulus 
less than 1, the general solution now takes the form 

(m) (m) 



U n+l U n+p 



with m = p — 2 and as before 



m I 



;=0 l<i 1 <i 2 <-..<ii<m t=l l<r<s<l 



Note the absolutely remarkable fact that we obtain again an expression in terms of the 
tau-function for the KP hierarchy, but with different kinematics, in the form of an implicit 
relation between the p's and <j's of eq. fl6.29Q . 

Imposing moreover the vanishing of the first terms U-i = U-i = ... = it_ p _|_i, we find 

that 

\i = x\ p ~ 1 " >m ~ 1 TT Si—Ei. ( 6 .32) 
A , A Pi - Pj 

This fixes completely the solution to our combinatorial problem, and in particular gives 
a compact expression for the generating function Rq for rooted bipartite p-valent planar 
graphs. For p = 4 for instance, it reads 

-SggiR 2 + 3g 2 g 2 R 4 
* = l-3 gjl g (633) 

while R satisfies R = 1 + 3ggiR s . 
6. 7. Constellations 

The solution for general p-constellations is similar to that for pure p-valent bipartite 
graphs. Indeed, we find that the general solution has the exact same form ( |6.30| ) as in the 
previous section, with given by ( |6.31j ) and pi,qi,Cij defined as in ( |6.29| ). The only 
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difference is that now m may take a larger value say (p — l)k — 1 in the case when only 
gr, gi, g 2l gk are non-zero (namely of constellations with white p valent vertices, and black 
vertices with valences p, 2p, ...,kp). The x's entering these formulas are now the generically 
distinct roots with modulus less than 1 of the linearized characteristic equation, itself a 
polynomial of degree (p — l)k — 1 of x + 1/x, reading 



x m-j(p-l) 



v 7 i>l m=0 j=0 v 7 v J 7 

(6.34) 

For illustration, in the case of 3-constellations with say only g and cj2 non-zero (R n 
obeys eq.( |5.12| ) with gi = 0), we get the linearized characteristic equation: 



1 - g 2 g 2 R 3 (x + - + 2)(x 2 + \ + x + - + 6) = (6.35) 

X X 1 X 



of degree 3 in x + 1/x. Picking the three roots with \xi\ < 1, we finally get the generating 
function for rooted 3-constellations with white trivalent and black hexavalent vertices 

1 - 17 g 2 g 2 R 3 + 25g 4 giR 6 
H ° = R l-WgJ (6 ' 36) 

while R satisfies the equation ( [4. 5| ) with only g and g 2 non-zero, namely R = 1 + lQgg 2 R' i - 
6.8. Ising model 

We have not been able to find a nice structure for the general solution of the even- valent 
bipartite case in general. In the particular case of the Ising model with zero magnetic field 
(eq. (|5.17[ )), we have been able to derive all possible solutions involving only one integration 
constant. 

The usual linearization of eq.( |5.17| ) yields the following characteristic equation 
1 c 1 -gV\ ( 1 c 1-gV 

X+- + —-TZ 7—7W ~ \X+ — — — ^— X 



x gV(l-3gV) gV J\ x gV(l-3gV) gV , o _. 

f 1 l + gV\ (6 ' 37) 
*{ X+ x + ^h)= 

of degree 3 in x + 1/x. As already discussed before, we would need in principle to find the 



full solution to (|5.17|) that converges at large n, including all three xi, x 2 , £3, respectively 



the roots of the three factors in eq. (|6.37| ) with modulus less than 1, and therefore including 
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also three integration constants. We now display the solutions with one x for each of the 
three factors in (|6.37|) . 



For both x = x\ and x = X2, we have found that 



while for x = x\\ 



and for x = X2'- 



V n = V UnUn+S , u n = l-Xx n (6.38) 

Un+lUn+2 

A(l - x)(l - x 2 )x n 
R n = R- V— - - — (6.39) 

Un+lUn+2 

R n = R + V X <- 1 - x)(1 - x2)xn (6.40) 

Un+lUn+2 



For x = X3 however, the solution is quite different: 
= y _ u n u n+3 



Un+lUn+2 

u n = l- 2Xx n - z\ 2 x 2n 



((c - 2)(x + J) + c - 8) ((c + 2){x + i) + c + 



(6.41) 



(a; 2 + ^-(c-4)(x + i)-(c-2))(x 2 + ^ + (c + 4)(x + i) + (c + 2)) 
The growing complexity of the solutions lets us expect a quite involved general three 
x-solution, yet to be found. 



7. Continuum limit 

The exact solutions of Sect. 6 allow us to investigate the properties of the correspond- 
ing classes of planar graphs in terms of the geodesic distance, in the limit when the latter 
becomes large. This limit must clearly be taken simultaneously with the so-called critical 
limit of large graphs, reached in turn by letting the various weights per vertex approach 
some critical locus, corresponding to approaching the finite radii of convergence of the 
various combinatorial series involved. 

The large n limit of R n = R+aiX™ + ... where x± is the largest root with modulus less 
than 1 of the characteristic equation, leads to the natural definition of a correlation length 
£ = — Log|xi|, governing the exponential decay of R n — R ~ exp(— n/£) as a function 
of the geodesic distance n. A good continuum limit may therefore be reached by letting 
x\ (and possibly other x's) tend to the value 1, while simultaneously keeping n/£, the 
continuum geodesic distance, fixed. This in turn implies certain relations between the 
vertex weights are reached, via the characteristic equation relating them to the x's. These 
relations express nothing but the abovementioned critical limit, which must therefore be 
taken simultaneously with the continuum one. 
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7.1. Tetravalent case 

We illustrate the above with the case of tetravalent graphs, with R n given by (|6.7|). 
The critical limit g — > g c = 1/12 is reached by taking say 

g = g c (l-e*), R = (7.1) 

1 + e z 



with R c = 2, and the solution of eq.( |6.2|) with modulus less than 1 reads 



l + 2e 2 -6y/3 (2 + 6 2 ) 
1 -e 2 



xiB) = /V^t^ (7 . 2) 



hence x = e ev ^ + 0(e 2 ) as e — > 0, and £ oc 1/e. Finally setting 

n = ^ (7.3) 

we may simply express the continuum limit of the quantity R n , or more interestingly that 
of R — R ni generating two- leg diagrams of tetravalent planar graphs with distance at least 
n between the two legs. We find that 

£ 6 K sinh 2 ( J\r 

This is the continuum two-point correlation function for random surfaces with two marked 
points at (rescaled) geodesic distance larger or equal to r. It coincides with the scaling 
function derived in jT7|], by use of transfer matrix formalism. This in turn yields the 



continuum two-point correlation function for random surfaces with two marked points at 
(rescaled) geodesic distance r: 



cosh I \/ %r 

G{r) = -F'(r) = 3^ VV ( (7.5) 

sinh 3 ( W | r 

This result may in turn be interpreted in terms of graphs with large but finite number 
iV of vertices. Indeed, the above scaling relations (|7.1| ) and (7^3) imply the following 



relation between the correlation length and the deviation from the critical point £ ~ 
(g c — g)~ v ', with the exponent v = 1/4, and therefore the fractal dimension dp = 1/v = 4 
for the present model of random planar surfaces. More concretely, this tells us in particular 
that the number of faces in a large graph lying at geodesic distance < n from the external 
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one behaves as n dp = n 4 . A simple measure of this number is indeed the ratio R n \ g N / Rq\ q n 
of the corresponding coefficients of g in the two power series expansions, giving the 
proportion of graphs with the two legs distant by at most n to that with the two legs in 
the same face. Using our exact solution and performing a saddle-point expansion, we find 

lim ^ikL „ A n 4 (7 6) 

n^oo Ro\gN 56 

giving an explicit illustration of the fractal dimension 4. 

This suggests to set n = rA 1 / 4 and to write R n \ g N, again by use of a saddle-point 
expansion, as 

4(12)* r °° 



R n\9» '--^3/2 



duu 2 (l + Rej zr (r v /Z m)) (7.7) 



and similarly for R\„n with T replaced by 0. The ratio R n \ g N / R\ g N gives the probability 
P(r) for a random surface with two marked points that their geodesic distance be less or 
equal to r: 

P(r) = —= \ duu 2 e- u2 (1 + ReF(r^iu) ) (7.8) 
Jo v 1 

7.2. Arbitrary even valences, multicriticality 

Repeating the analysis of previous section for the case of graphs with even valences say 
up to 2m + 2, various critical points may be reached in the space of weights gj. Concretely, 
picking the particular weights that ensure the following form for (|2.3| ): 



9c ~9 V c -V 



9c V V c 



(7.9) 



where V = gR, g& = g, g%j = gi Zj fixed by the form ( |7.9[ ) and V c = to/6, g c = m/ (6(to+1)), 
the corresponding multicritical limit is obtained by setting 

9 = 9c(l-e^), ^ = i?c 1 _ 1 ~f +1) (7.10) 
with R c = V c /g c = to + 1. Remarkably, the characteristic equation ( |6.9| ) then turns into 



Xm(x) = ( \^-) m P m 0-/- (x + \ - 2) ) = (7.11) 



for some fixed degree m polynomial 

m 

(2Z + 1)! (m-l)l 



m j. . 

p m (u) = £("">' r„^, fm _m < 7 - 12 ) 



Z=0 
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This allows to get the leading behavior of the various x's as 

Xi = e -^ + 0(e 2 ) (7.13) 

where a 2 are the roots of P m , and aj are taken with positive real part. Consequently all 
the x's tend to 1 simultaneously, and the correlation length of the problem reads £ ~ 1/e, 
hence we now have the relation £ ~ (g c — g)~ v with v = 2(m + 1), i.e. a fractal dimension 
d F = 2(m + 1). 

The multicritical continuum limit is therefore still obtained by setting ( [7.3| ), and the 
exact solution ( |6.10| ) yields the following two-point correlation of random surfaces with 
multicritical weights (known to simulate non-unitary matter conformal field theories with 
central charges c(2, 2m + l) = 1 — 3(2m — l) 2 /(2m + l) coupled to two-dimensional quantum 
gravity), with two marked points at geodesic distance less or equal to r: 

d? TV T 

^( r ) = -2^LogW(sinh(a 1 ^ ,sinh(a 2 ^ , sinh (a m 0) (7.14) 

where W(/i, f2, f m ) stands for the Wronskian determinant det [f^ 1<i - <m - 

It is known from matrix model solutions that the general case of arbitrary valences 
leads to the same multicritical points. In particular we expect the scaling functions ( |7.14j ) 
to be the same at these points. The same remark applies to constellations as well. New 
multicritical points corresponding to conformal theories with central charges c(p, q) = 
1 — 6(p — q) 2 /(pq) for p, q two coprime integers can be reached within the framework of 
two-matrix models, corresponding to the general bipartite graphs. In the latter case, we 
expect some new scaling functions, characteristic of these other universality classes. An 
example will be given in next section, when discussing the Ising model. 

7.3. Critical/continuum limit in general 

The (multi-) critical continuum limits are reached by letting a number of the roots x 
of the characteristic equation tend to 1 simultaneously. 

An alternative route for deriving the critical continuum limit of say the tetravalent 
case would have been to postulate the form R n = R(l — e 2 J r (ne)) and plug this ansatz 
into the recursion relation (|5.1|). With g = g c {l — e 4 ), expanding the equation up to order 
4 in e, we arrive at the following differential equation for T: 

T"{r) -3^ 2 (r) - 6F(r) =0 (7.15) 
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The function T of eq.( [7.4| ) is the unique solution to (|7.15| ) such that TifS) = oo and 
.F(oo) = 0. 

More generally, we may derive differential equations for the multicritical cases of 
Sect. 7. 2 as well. These take the form 

K m+1 [l + F] =n m+1 [l] (7.16) 

where 7£ m [/u] is the mth KdV residue (d 2 — w,) m-1 / 2 |d-i where d = d/dr [EDI. For instance, 
in the case of multicritical tetra/hexavalent graphs (m = 2) eq. (|7.16|) reads 

^^^-lOTi^T'^^-lO^'^^-bi^'ir^ + lOi^ir^ + ^Oi^ir^ + ^O^ir) = (7.17) 

Even more generally, recall that the KdV residues naturally arise (e.g. in the context 
of matrix model solutions) when solving the differential operator equation [P, Q] = 1, where 
Q = d 2 — u and P some degree 2m + 1 differential operator. The equation indeed boils 
down to 2d/dr(TZ rn+ i[u\) = 1. In the present case, we rather have to write the equation 
[P, Q] = 0, which turns into 7?. m _|_i [it] =const. 

At the discrete level, comparing say (|5.1| ) with the equations determining the matrix 
model solution for tetravalent graphs of arbitrary genus, we simply would have to replace 
1 by n/N in the r.h.s. of (|5.1| ). We may think of the differential operator Q = d 2 — u as 
the continuum limit of the operator Q (|5.2|) defined in Sect. 5. 2, now acting on functions 
of the variable r = ne. In the other cases described in Sect. 6, we always have such an 
operator Q at hand and again the recursion relations resemble strongly those obtained 
in the solutions of the multi-matrix models describing two-dimensional quantum gravity 
coupled with matter, up to the same substitution 1 — > n/N. This suggests that the relevant 
(coupled) differential equations governing the (multi-) critical continuum limit should read 
[P, Q] = 0, with Q the continuum limit of the operator say Qi used in our approach, 
taking the form of a differential operator of degree q say, and P a differential operator 
of degree p coprime with q. Our claim is that the generalized two-point functions .F(V) 
for random surfaces in the presence of critical matter (corresponding to confer mal field 
theories with central charges c(p, q) < 1), with two marked points at geodesic distance less 
than r, should be governed by [P, Q] = 0, where Q = d q — qud q ~ 2 .., and u = 1 + T. This 
is illustrated in the case of the Ising model in the next section. 
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7.4- Ising model 

The multicritical limit of the Ising model is obtained as follows. Starting from the 
equations ( fO| ) rewritten as W(R) = 0, W a polynomial of degree 5, the tri-critical points 
are determined by setting W'(R) — W"(R) = 0, and we find that c c = ±4 while g c = 10/9, 
and R c = —3/5, V c = —3/10. The tricritical limit is approached by setting 

16 

c = 4, R = R c (l-e 2 ) g = g c (l-—e 6 ) (7.18) 

We note that in the characteristic equation ( |6.37| ) only the first and last factor tend to 
x + 1/x — 2 as e — > 0, while the middle one tends tox + l/x + 10. This means that only 
two of the three x's tend to 1 in this limit. More precisely, we have xs = e~^ e + 0(e 2 ) 
and X\ = e _2v ^ e + 0(e 2 ). This displays the fractal dimension of graphs with critical 
Ising configurations, namely dp = 6, obtained by expressing the correlation length £ ~ 
(g c — g)~ v ', where v = 1/6. Recall however that the distance n or its rescaled version r 
are slightly different from the true geodesic distance in the Ising tetravalent graphs, and 
the fractal dimension measured here might be different from that associated to the true 
geodesic distance. 

Further substituting 

R n = R(l - e 2 J r (ne)), r = ne 

R n (7-19) 



V 

v n 



c + g(R n +i + R n + R 



n— 1, 



into the recursion relations ( |5.17|) together with Q7.18 ), we find by expanding up to order 



6 in e that the two-point function T obeys the differential equation 

j^ 4 ) (r) - 18^)^" (r)- 18^" ^ = (7.20) 

This equation is precisely what one would get by writing [P, Q] — for differential operators 
Q = d 3 - 3ud - 3u'/2 and P of order 4, u = 1 + T. 

Looking for convergent solutions in the form T{r) = ae~ kr ', we find that k\ = or 
k<2 = 2\/3, corresponding to e~ kir = x% and e~ k2l ~ = respectively. Proceeding like in 
the discrete case, we may now solve the differential equation order by order in e~ kir , with 
a double power series expansion •F(r) = p>0 a mjP e~ r (' mkl+pk2 \ in terms of the two 
integration constants A = a^o, u = ao,i, while ao,o = 0. The differential equation ( |7.20j ) 
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indeed simply amounts to a recursion relation on the coefficients a m>p . Resumming the 
series for J-{r) finally yields: 

JF( r ) = --^-Log ( 1 - -e~ rVE - A e -2rV3 _ >?_„-2rVS 



dr 2 & V 6 12 288 

72 VC + 3456 A ^ J 

The two integration constants are further fixed by the boundary conditions. The latter are 
obtained by requiring that the recursion relations (|5.17| ) also make sense at n = 0, 1, namely 
R-i = V-i = 0, while lim n ^o V n -iV n -2 = 0. The first conditions give .F(O) = oo, while 
the latter implies R-1R-2 = 0. As in the case of planar graphs with only tetra/hexavalent 
vertices, this implies a higher order vanishing of the argument of the logarithm in ( |7.21|) , 
which plays the role of continuum limit of u n , while the condition implies that both u-\ 
and U-2 vanish. This is easily solved into A = —12(17 + 12\/2) and fi = 12(4 + 3\/2), so 
that finally 

d? / 

;F(r) = --^Logfsinh (r(V6+V3)) +(17+12^) sinh (r(VE-VS)) -2(4+3^2) sinh(rV3) 

(7.22) 

and we also get the correlation Q(r) = — jF'(r) for the Ising model on random surfaces, 
with two marked points at (special) geodesic distance r. 



8. Conclusion 

In this note, we have addressed the problem of enumeration of various types of pla- 
nar graphs with external legs, while keeping track of a suitable geodesic distance between 
these legs. The basic tool we used are blossom-trees, namely trees carrying the minimal 
information needed to close them back into planar graphs. This information is essentially 
contained in the two types of leaves (black and white), and we have devised a compact 
algebraic way of keeping track of the geodesic distance between legs, by introducing opera- 
tors Q describing the structure of the rooted trees around their first vertex. This operator 
acts formally on a basis |n) indexed by relative integers, and may as well be viewed as 
acting on sequences {p n }nezi y i a the shift operator a and its relative integer powers, and 
a number of (combinatorial) diagonal operators. This definition is clearly borrowed from 
that of the Q operator of matrix models, that generate the multiplication by an eigenvalue 
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A on the basis of (right) monic (bi-) orthogonal polynomials p n (X), n = 0,1,2... In this 
framework, the main recursion relations for the coefficients of Q are obtained by consider- 
ing the operator P, acting on the p n (A) by differentiation w.r.t. A. The canonical relation 
[P, Q] = 1 determines in fact all the functions of the problem, and turns into differential 
equations in the critical scaling limit. Here, by analogy with this case, we have been led to 
set [P, Q] = 0, for some operator P still awaiting a good combinatorial meaning. Neverthe- 
less, this equation also turns into differential equations for the physical quantities of our 
problem. The mysterious part of this correspondence is that the natural variable in the 
matrix model approach is the rescaled cosmological constant, which generates the topolog- 
ical expansion of the free energy, namely the expansion in powers of r of the function u(r) 
has coefficients corresponding to planar graphs of fixed genus. This must be contrasted 
with the present findings, where r has the meaning of geodesic distance. This seems to 
indicate that a more general structure should exist, that includes both the topological and 
geodesic directions, probably some suitably defined matrix model of sorts. 

In view of these strong analogies with matrix model solutions, we may want to char- 
acterize the class of possibly decorated planar graphs for which a tree formulation exists as 
that for which there exists a matrix model formulation admitting a solution via orthogonal 
polynomials. This would exclude for instance the case of the three-state Potts model, a 
generalization of the Ising model whose configurations are graphs with vertices of three 
possible colors, and edge weights w a ^ = e KSa ' h according to the colors a, b of the adja- 
cent vertices. More generally, loop models on graphs also correspond to matrix models 
without orthogonal polynomial solutions: their configurations are simply mutually- and 
self-avoiding loops drawn on the edges of planar graphs, with a weight n per loop (the 
so-called O(n) model, extensively solved in ^TJ). It would be extremely interesting to 
investigate any of these models using tree techniques. 

A final striking outcome of our work is the emergence of soliton-like tau functions en- 
tering the explicit exact formulas for a number of generating functions for two-leg diagrams 
with legs distant by at most n. A reason why this should happen in the first place may 
be perhaps traced back to the integrability of the recursion relations we have obtained. 
A first indication of this integrability is the existence of "integrals of motion" for these 
equations. For instance, the recursion relation of the tetravalent case, eq.([5J]), has the 
following integral of motion: 

f(R n ,R n+1 ) = const. 

(8.1) 

f(x, y) = xy(l - gx - gy) - x - y 
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as is immediately checked by forming 



f(R n , Rn+l) — f{Rn, Rn-l) ~ (Rn+1 ~ Rn-l){Rn — 1 — gR n (R n +l + R n + Rn-l)) (8-2) 

More generally, one may construct such integrals of motion for all the models studied in this 
note. Again, it may be that this integrability property relates to the existence, for the same 
planar graphs but with arbitrary topology rather than planar with fixed geodesic distances, 
of matrix model formulations that are solvable via orthogonal polynomial techniques. 
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